Magnetic field generation from non-equilibrium phase transitions. 



D. Boyanovsky* 

Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA 
and LPTHE, Universite Pierre et Marie Curie (Paris VI) et Denis Diderot (Paris VII), 
Tour 16, ler. etage, 4, Place Jussieu, 75252 Paris, Cedex 05, France 

H. J. de VegaQ 

LPTHE, Universite Pierre et Marie Curie (Paris VI) et Denis Diderot (Paris VII), 
Tour 16, ler. etage, 4, Place Jussieu, 75252 Paris, Cedex 05, France 
and Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA 

M. SimionatoQ 

£\j , Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA 

Q ■ (Dated: February 1, 2008) 

■ We study the generation of magnetic fields during the stage of particle production resulting from 
spinodal instabilities during phase transitions out of equilibrium. The main premise is that long- 
wavelength instabilities that drive the phase transition lead to strong non-equilibrium charge and 
current fluctuations which generate electromagnetic fields. We present a formulation based on the 

j /^ ' non-equilibrium Schwinger-Dyson equations that leads to an exact expression for the spectrum 

of electromagnetic fields valid for general theories and cosmological backgrounds and whose main 
ingredient is the transverse photon polarization out of equilibrium. This formulation includes the 
dissipative effects of the conductivity in the medium. As a prelude to cosmology we study magneto- 
genesis in Minkowski space-time in a theory of N charged scalar fields to lowest order in the gauge 
coupling and to leading order in the large N within two scenarios of cosmological relevance. The 
long-wavelength power spectrum for electric and magnetic fields at the end of the phase transition 
is obtained explicitly. It follows that equipartition between electric and magnetic fields does not 
hold out of equilibrium. In the case of a transition from a high temperature phase, the conductiv- 
qq ■ ity of the medium severely hinders the generation of magnetic fields, however the magnetic fields 

generated are correlated on scales of the order of the domain size, which is much larger than the 
, magnetic diffusion length. Implications of the results to cosmological phase transitions driven by 

■ spinodal unstabilities are discussed. 
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I. INTRODUCTION 

There is compelling astrophysical evidence for the existence of magnetic fields in the Universe . Recent advances 
in observational techniques mainly through Faraday rotation (RM) complemented with an independent measurement 
of the electron column density via pulsar dispersion measurements (DM) or X-ray emission from clusters]^] indicate 
that the strength of these astrophysical magnetic fields is of order of fiG and are coherent on scales up to that of galaxy 
clusters 100 Kpc — 1 Mpc|l], ^]. These magnetic fields are now deemed to play an important role in the evolution and 
dynamics of galaxies but their origin is still largely unknown ||, ||, ||, 

The galactic dynamos are some of the promising mechanisms to amplify pre-existing seeds ^ [?], In its simplest 
conception a linear or kinetic dynamo transfers energy from differential rotation to the build-up of a magnetic field with 
a typical growth rate determined by the rotation period of a protogalaxy ~ Gyr~ 1 . There are many alternative versions 
of dynamo theory, and some of the most promising require turbulent flows Dynamos amplify seeds but an initial 

seed must be present. The proposals to explain the initial seeds can be classified as astrophysical or of primordial 
cosmological origin. An important astrophysical mechanism is the Biermann battery (for a recent discussion see |J) 
which relies on gradients in the electron number density and pressure which arc in different directions that can arise 
behind hydrodynamic shocks^]. The primordial cosmological mechanisms purport the generation of seed magnetic 
fields during different stages of the Early Universe^). As originally observed by Turner and Widrow[^[ |J the coupling 
of the electromagnetic field to gravity is conformally invariant, resulting in that cosmological expansion per se will not 
generate primordial magnetic fields, the coupling of charged fields to gravity is in general not conformally invariant 
and can lead to magnetogenesis through the minimal coupling to electromagnetic fields. These authors studied a wide 
range of possibilities for primordial magnetogenesis with encouraging results. More recently a thorough study of the 
generation of (hyper) magnetic fields during inflation concluded that the amplitude of the seeds on cosmologically 
relevant scales is probably too small[jl(|. These authors included dissipative effects in the medium via a kinetic 
approach that includes the conductivity. 

The amplification of electromagnetic fluctuations during inflation, by the end of inflation and from inflation to 
reheating has been studied by several authors in the context of different inflationary models [0 (see criticisms in 
ref. jj3). Dolgov [jl3| suggested that the breaking of conformal invariance through the trace anomaly may also lead to 
significant magnetogenesis (see also Jl4|]) and non minimal electromagnetic couplings has also been considered || 0, |l5| . 
Other proposals suggest that primordial magnetic fields can be generated much later during the electroweak phase 
transition! 16| or in connection with dark energy during the epoch of large scale structure jlTj . For a pioneering proposal 
for magnetogenesis from cosmological first order phase transitions see ref. [[l8j . 

While certainly there is no dearth of proposals an important aspect that is clear from the present body of work is 
that a consistent framework to study the generation of primordial magnetic fields must include consistently (or semi 
phcnomenologically) the dissipative effects associated with the high conductivity of the medium. This aspect was 
already highlighted in the seminal work of Turner and WidrowQ and again by Giovannini and Shaposhnikov|l(| ^2| . 

This is clearly more relevant if the medium is a hot and dense plasma as prevailed during reheating, preheating 
and most likely any phase transition associated with particle physics scales. 

The goals of this article: In this article we begin a program to assess the generation of magnetic fields through 
the non-equilibrium processes associated with supercooled (second order) phase transitions. Here we focus on the 
process of spinodal decomposition, namely the process of phase separation, domain formation and growth resulting 
from a non-equilibrium phase transition as the mechanism that produces strong charge fluctuations which lead to 
magnetogenesis. Although we focus primarily on spinodal decomposition, similar arguments will apply for very weak 
first order phase transitions with small latent heat and nucleation barriers, since in this case nucleation will be almost 
indistinguishable from spinodal decomposition. We study models in which there is a charged scalar sector coupled to 
an abelian (hypercharge) gauge field without breaking the U(l) symmetry associated with the (hyper) charge. The 
charged scalar and abelian gauge fields are supposed to be part of a larger multiplet of fields pertaining to a larger 
underlying theory. 
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We develop the framework to calculate the power spectrum of magnetic fields generated out of equilibrium directly 
from the underlying quantum field theory. An exact expression for the spectrum of magnetic fields generated through 
non-equilibrium processes is ob taine d directly from the non-equilibrium Schwinger-Dyson equations for the transverse 



gauge field propagator [see eq.( 3.45 )1- This expression is general and valid for all types of charged matter fields, the 
main ingredient being the non-equilibrium transverse polarization for the (hypercharge) gauge field. 

This framework allows to include the effects of a conductivity in the medium and is a generalization of a previous 
study of photon production out of equilibrium Jl9[ . 

After establishing a connection between the dynamics in cosmological space times and the simplified case of flat 
Minkowski space time with a model for the phase transition, we apply these methods to study the power spectrum 
of magnetic fields during quenched phase transitions. In order to study the non-perturbative aspects of spinodal 
decomposition, such as the growth of unstable modes below the phase transition we invoke a model of N charged 
fields in the large N limit. Most particle physics models (GUT's, SUSY, etc.) contain a large number of scalar fields, 
thus warranting the use of the large N approach. 

The purpose of this study is to extract the important and robust features that lead to the generation of magnetic 
fields from these non-equilibrium processes in a simpler setting and to build intuition into the main physical aspects. 

We postpone the application of these methods to the full cosmological situation to forthcoming articles. 

As mentioned above, many mechanisms for primordial magnetogenesis have been studied previously Our 
main goal here is to provide reliable quantitative estimates for primordial magnetogenesis during out of equilibrium 
phase transitions. In particular we focus on the dynamics of magnetogenesis in a minimally coupled scalar-gauge theory 
where the non-equilibrium dynamics is a consequence of spinodal instabilities. The formulation presented here is very 
general, it includes consistently absorptive effects such as the conductivity and can be implemented in other scenarios. 

In section II we present the main physical picture behind the mechanism of magnetogenesis studied. In section III we 
introduce the model under consideration and discuss the issue of gauge invariance. Section IV discusses the similarity 
between Friedmann-Robertson- Walker cosmologies and flat Minkowsky space time that warrants a preliminary study 
of magnetogenesis within the simpler setting. As a prelude to cosmology, section V briefly reviews the main aspects of 
non-equilibrium dynamics of spinodal decomposition in the scalar sector in Minkowsky space time, which are present 
in the cosmological setting. In section VI we obtain the exact expression for the power spectrum of magnetic fields 
directly from the non-equilibrium Schwinger-Dyson equations for the transverse photon propagator. We compare the 
non-equilibrium result to the more familiar equilibrium setting. 

In section VII we implement this framework to study magnetogenesis in two cosmologically relevant settings that 
purport to model a phase transition during inflation and a radiation dominated era. We discuss in detail how the 
conductivity of the plasma emerges in this formulation and obtain the energy density on long-wavelength magnetic 
fields extracting features that will remain in the full cosmological problem. 

II. THE PHYSICAL PICTURE 

Whether a symmetry breaking phase transition occurs in local thermodynamic equilibrium (LTE) or not is a detailed 
question of time scales. Two important time scales must be compared, the time scale of cooling through the critical 
temperature r coo / = T(t)/T(t) and the relaxation time scale of a fluctuation of typical wavelength k 7 r(k). If cooling 
through T c occurs on time scales slower than the relaxation scale, i.e, t coo i 3> t(&) then a fluctuation of wavevector k 
adjusts to LTE and for these scales the phase transition occurs in LTE. On the other hand if r coo / <C r(k) fluctuations 
of wavelength cannot adjust to the conditions of LTE and fall out of equilibrium, i.e, freeze out and for these 
scales the phase transition occurs on very short time scales like a quench from a high temperature phase into a low 
temperature phase p9 . While short wavelength modes typically have very short relaxation time scales and remain in 
LTE during the phase transition, long-wavelength modes undergo critical slowing down near T c and their relaxational 



dynamics becomes very slow pi 
wavelength modes freeze out fal 



. Thus as the temperature nears the critical with a non-vanishing cooling rate, long- 
ing out of LTE, namely for these modes the phase transition occurs out of equilibrium. 
The non-equilibrium dynamics below the critical temperature is described by the process of spinodal decomposition: 
long- wavelength fluctuations become unstable and grow initially exponentially in time (in Minkowski space-time) . The 
field becomes correlated within regions characterized by a time-dependent correlation length £(i) and the amplitude 
of the long-wavelength fluctuations becomes non-perturbatively large since the mean-square root fluctuation of the 
field will begin to probe the broken symmetry states with (cf)' 2 ) oc 1/A with A the scalar self-coupling. This is the 
hallmark of the process of phase separation, the correlated regions correspond to domains, inside these domains the 
field is near one of the vacuum states. A particular time scale, the nonlinear time t n i oc ln(l/A) , determines the 
transition from a regime of linear instabilities to one in which the full non-linearities become important. This time 
scale roughly corresponds to when the mean square root fluctuation of the scalar field samples the broken symmetry 
states and the phase transition is almost complete. At this stage the amplitude of the long-wavelength fluctuations 
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become of 0(1/A), the phases freeze out and the long-wavelength fluctuations become classical but stochastic p2[-p3[. 

Consider now the situation where the scalar field carries an abelian charge and is coupled to the electromagnetic 
gauge field. The strong spinodal fluctuations in the scalar field will induce fluctuations in the current-current correla- 
tion function, and while the expectation value of the current must vanish by translational and rotational invariance, 
the current and charge correlators will have strong non-equilibrium fluctuations. These current fluctuations will in 
turn generate a magnetic field with a typical wavelength corresponding to the wavelength of the spinodally unstable 
modes. 

This is the main premise of this article: the spinodal instabilities which are the hallmark of a non-equilibrium 
symmetry breaking phase transition will lead to strong charge and current fluctuations of the charged scalar fields 
which in turn, lead to the generation of magnetic fields through the non-equilibrium evolution. 

Thus while the mechanism that drives magnetogenesis is the same as in most other scenarios, namely, the ampli- 
fication of charge fluctuations JToj, |l7| , in this scenario, the amplification is the result of the spinodal instabilities of 
long-wavelength fluctuations during the phase transition. 

The formation and growth of long-wavelength correlated domains entails strong non-equilibrium charge and density 
fluctuations. These would not be possible if there were long-range interactions, however in a high temperature plasma 
the long-range Coulomb interaction is actually screened. In a high temperature plasma the Debye screening length 
is £d oc (eT)" 1 !^. Thus long- wavelength fluctuations will not be hindered by the Coulomb interaction which is 
effectively screened over very short distances in a high temperature plasma. 

Primordial magnetic field generation by an electric charge asymmetry was previously studied by Dolgov and Silk |25f| 
assuming an early stage in which electromagnetism was spontaneously broken. In that scenario long-range forces are 
screened by the effective mass of the gauge fields arising from the expectation value of the Higgs field. At high 
temperature there is also screening from the plasma, which, however, was not taken into account inp5|. 

The situation that we study is rather different, the gauge symmetry is never spontaneously broken, but long range 
forces are screened on scales given by the Debye screening length £d ~ (eT) -1 . 

The main ingredients that must be developed in order to understand the generation of magnetic fields through this 
non-equilibrium process are: 

• A consistent framework to compute the spectrum of generated magnetic field, namely (B(k,t) ■ B{—k,t))/V 
with B(k,t) the spatial Fourier transform of the Heisenberg magnetic field operator and V the volume of the 
system. 

• We anticipate that plasma effects must be included to assess the generation and eventual decay of the magnetic 
fields. If there is a large conductivity in the medium the magnetic field will diffuse but also its generation will 
be hindered. This will be a point of particular importance within the cosmological setting j^, flo[ . 

These ingredients will be analyzed in detail below. 

III. THE MODEL 

The cosmological setting, in which we are primarily interested, corresponds to either a symmetry breaking phase 
transition during (or leading to) an inflationary stage or in a radiation dominated Universe. Such phase transition 
is in principle different from the electroweak one 1 and presumably occurs at a much higher energy scale, such as 
the GUT scale ~ 10 15 Gev but is assumed to be described by a particle physics model that includes many fields with 
(hyper)-charge either fermionic or bosonic. We will not attempt to study a particular gauge theory phenomenologically 
motivated by some GUT scenario, but will focus our study on a generic scalar field model in which the scalar fields 
carry an abelian charge, which for simplicity we will take to be electric charge. 

We first focus on the generation of magnetic fields through the non-equilibrium dynamics of the phase transition in 
Minkowski space-time, and we will argue in a later section that the case of Friedmann-Robertson- Walker cosmologies 
in conformal time is very similar to the Minkowski space-time case. Therefore most of the technical aspects developed 
within the framework of Minkowski space-time will be translated vis-a-vis to FRW cosmologies. 

A related study was carried out in reference Jigf] where dynamical aspects of screening and photon production were 
studied both in the case of non-equilibrium symmetry breaking as well as of parametric resonance. 



If the electroweak phase transition is weakly first order, nucleation will be almost indistinguishable from spinodal decomposition and 
the phenomena studied here may be of relevance. 
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While some aspects of the dynamics of magnetic field generation are related to those studied in |l£j], there are 
important differences that warrant a detailed treatment. In particular, the results of photon production out of 
equilibrium could be used if equipartition held, since both electric and magnetic fields contribute to photon production. 
However, as it will be studied in detail below, out of equilibrium equipartition is no longer valid and the generation 
of magnetic fields must be studied in general, including the effects of the conductivity. 

While details may depend on the specific models, we seek to extract the robust features of the phenomena which 
will be captured by a model of scalar electrodynamics with N charged scalar fields. [Within the framework of unified 
theories, this abelian gauge field should be interpreted as carrying a hypercharge quantum number]. Furthermore, the 
scalar sector contains also a neutral field which plays either the role of the inflaton or some generic field that rolls to 
the minimum of the potential leading to symmetry breaking, in such a way that electromagnetism is not spontaneously 
broken to describe the correct low energy sector with unbroken U(1)em- We will take the neutral and the N complex 
(charged) fields to form a scalar multiplet under an 0(2N + 1) isospin symmetry. As the neutral field acquires an 
expectation value this isospin symmetry is spontaneously broken to 0(2N). Since by construction only the neutral 
field acquires a non vanishing expectation value under the isospin symmetry breaking the photon remains massless 
(it will obtain a Debye screening mass from medium effects). 

The model under consideration is defined by the following Lagrangian density, 

C = \d^a d»a + D^D^ - \f^F^ - m 2 Qa 2 + ^ _ A fi^ + ^ 
where to simplify the expressions we have introduced the following notation 

N N 
r—1 r—1 

Furthermore, anticipating a non-perturbative treatment of the non-equilibrium dynamics of the scalar sector in a 
large N expansion, we have rescaled the quartic coupling in such a way as to display the contributions in terms of 
powers of l/N, keeping A fixed in the large N limit. 

Gauge invariance: Gauge invariance must be treated very carefully in an out of equilibrium situation. In 
this abelian theory the constraints arising from gauge invariance can be accounted for in a systematic and fairly 
straightforward manner leading to a gauge invariant description as follows (for details see p6|). There are here two 
first class constraints 

n «0 ; V-E-p^O (3.1) 

where IIo is the canonical momentum conjugate to the time component of the vector field and p is the charge density. 
The second one is identified with Gauss' law, the generator of time independent gauge transformations. At the 
quantum level, these constraints annihilate the physical states simultaneously. In a Schrodinger representation of the 
quantum field theory, the first constraint implies that the wave-functionals do not depend on A$ . The condition that 
Gauss' law annihilates the physical states implies that the wave-functionals are only functionals of particular gauge 
invariant combination of the fields described in detail in ref . p(l| . In a Hamiltonian description this procedure is similar 
to the Coulomb gauge but we haste to emphasize that this procedure is fully gauge invariant. The instantaneous 
Coulomb interaction can be replaced, at the level of the path integral, by an auxiliary field Aq which is the Lagrange 
multiplier that enforces Gauss' law and should not be identified or confused with the temporal component of the 
vector field. The final result is that this gauge invariant description is equivalent to using the following Lagrangian 
density (for details see pp| ) 

C = \d»o d»a + e^t + lg^A T ■ d»A T + i(VA ) 2 - m 2 Qa 2 + - A + $t $ 

-iel T ■ ($ f V$ - V$ f $) - e 2 (A T - A%) - ie A ($$ f - (3.2) 

where $ is a gauge invariant local field which is non-locally related to the original fields, and At is the transverse 
component of the vector field (V • At = 0) and A$ is a non-propagating field as befits a Lagrange multiplier, its 
dynamics is completely determined by that of the charge density [|26| . 

The main point of this discussion is that the framework to obtain the power spectrum of the generated magnetic 
field presented below is fully gauge invariant. 
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IV. MAGNETIC FIELDS IN FRIEDMANN-ROBERTSON- WALKER COSMOLOGY 

The action for a charged self-interacting scalar field 4> coupled to an abelian (charge or hypercharge) gauge field A 
in a general relativistic background geometry is given by 



d X\J — g 



,2 it 



— ( 

AN 



1 



where 



A Fricdmann-Robcrtson- Walker line clement 



(4.1) 



(4.2) 



ds 1 



dt 2 - a 2 (t) dx 2 



(4.3) 



is conformally related to a Minkowski line element by introducing the conformal time r\ and scale factor C(rj) as 



dr\ 



dt 



C( V ) = a(t( V )) 



In terms of these the line element and metric are given by 

ds 2 = C 2 ( V ) (d V 2 -dx 2 ) ; 9liv = C 2 (r)) Vmv 
where if"" = diag(l, —1, —1, —1) is the Minkowski metric. Introducing the conformal fields 

A l ,(r),x) = A(t(r)),x) , $(r},x) = C{rj) <t>(t(r}),x) 
and in terms of the conformal time, the action now reads 



(4.4) 



(4.5) 



5' 



dr] d x 



rj^D^Dv® - M 2 {rj) 



A 

AN 



($*$) 2 



Ffiu F a0 rf v ^ 



with 



M 2 (ii) = m 2 C 2 (r]) 



C"(y) 
C(rj) 



(4.6) 

(4.7) 
(4.8) 



and the primes refer to derivatives with respect to conformal time. Obviously the conformal rescaling of the metric and 
fields turned the action into that of a charged scalar field interacting with a gauge field in flat Minkowski space-time 
with the scalar field acquiring a time dependent mass term. Here we have omitted the contribution from the conformal 
anomaly, which has been studied in ref . [fl3f . In particular, in the absence of the charged scalar field and neglecting 
the contribution from the conformal anomaly, the equations of motion for the gauge field A are those of a free field in 
Minkowski space time. This is the statement that gauge fields are conformally coupled to gravity and no generation 
of electromagnetic fields can occur from gravitational expansion alone without coupling to other fields or breaking the 
conformal invariance of the gauge sector. The generation of electromagnetic fields must arise from a coupling to other 
fields that are not conformally coupled to gravity (like the scalar field), or by non minimal electromagnetic couplings 
that would break the conformal invariance of the gauge fields. 

The conformal electromagnetic fields £, B are related to the physical E, B fields by the conformal rescaling 



E 



c 2 {n) 



B = 



B 



C 2 (v) 



(4.9) 



corresponding to fields of scaling dimension two. 

Since the action (4.6) is essentially that of coupled charged scalar fields with gauge fields in Minkowski space time 
(with a time dependent mass term for the scalars), much will be learned about the cosmological problem by first 
studying the main dynamical aspects in Minkowski space-time. 
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V. PHASE TRANSITIONS IN MINKOWSKI SPACE TIME: A PRELUDE TO COSMOLOGY 

Our ultimate goal is to study the generation of primordial magnetic fields during phase transitions either in an 
inflationary or radiation dominated eras which are the most likely scenarii for particle physics phase transitions. 

During inflation, the phase transition that we consider corresponds to the dynamics of small or large field models 
where the scalar field has a symmetry breaking or an unbroken potential, respectively. The rolling of the field 
corresponds to the expectation value evolving in time towards the minimum of the potential but also quantum 
fluctuations growing through the spinodal or parametric instabilities p2| , p7| , |28|. The initial state is generally 
described by a gaussian wave-functional j|7], ^] . 

In a radiation dominated cosmology, the initial state is that of local thermodynamic equilibrium at an initial 
temperature T » T c . The cosmological expansion leads to cooling T(t) — T/a{t) and at some given time the 
temperature equals the critical and the phase transition occurs. Using finite temperature field theory in an expanding 
background geometry, it is shown that the effective time dependent mass term includes the time dependent temperature 



corrections (for details see for example 
While in a cosmological background t 



29j and references therein). 

;he temperature falls off with the scale factor and the transition from the high 
temperature phase to the low temperature phase is driven by the expansion, in Minkowski space time we must model 
the transition. As described above, while short wavelength fluctuations are in LTE, long wavelength fluctuations will 
undergo critical slowing down and for these the phase transition will occur as a temperature quench from the high to 
the low temperature phase. 

A rapid (quenched or supercooled) symmetry breaking phase transition can be modeled in Minkowski space-time 
by a time-dependent mass term which changes sign suddenly from positive to negative |2^, |2^| . This quench approxi- 
mation can be relaxed by allowing a mass term that varies continuously with time with results that are qualitatively 
similar |3(| to the quench scenario. 

We consider two different scenarios that we invoke to model a quenched phase transition either during an inflationary 
era or a post-inflation radiation dominated era: 

i) : A vacuum phase transition, in which the initial state corresponds to a vacuum state of a free scalar theory with 
a positive squared mass term m 2 > 0. This corresponds to an initial pure Gaussian density matrix describing a free 
field theory of scalars with a mass m, L p^, 5s|, this initial state is evolved in time with the Hamiltonian with the 
symmetry breaking potential. 

ii) : A transition during a radiation dominated era. In this scenario, before the phase transition (t > to) the system 
is in the unbroken high temperature phase and the scalar fields acquire a positive thermal mass = (A/24 + e 2 /8)T 2 
arising from scalar tadpole and the gauge contributions to the scalar self-energy, whereas after the transition (t 3> to) 
the system is in the low temperature phase with negative mass square m 2 — —^i 2 . 

Thus we model the time-dependence of the mass term as 

m 2 (t) = M 2 9(-t) - fi 2 6{t) . (5.1) 
where we have chosen the transition to occur at t = and with M 2 > and given by 



^ + tor the radiation dominated case 

Detailed analytical and numerical studies of the dynamics of phase transitions in expanding cosmologies reveal that the 
qualitative features of phase separation and spinodal decomposition are well captured by this simple approximation |2^ , 
p?] . We postpone to a forthcoming article the full study of primordial magnetic field generation in a cosmological 

background. 

We emphasize that since the U(l) abelian vector field A^, namely the photon, couples to the charged scalar fields 
and the 0(2N) symmetry of the charged sector is unbroken, there is no Higgs mechanism and the photon does not 
acquire a mass. 

A. Scalar fields dynamics 

For completeness and to highlight the aspects of the non-equilibrium dynamics most relevant to the generation of 
magnetic fields, we summarize the main features of scalar field dynamics. For further details the reader is referred 

As described above the non-equilibrium evolution of long-wavelength modes begins with the spinodal instabilities 
which in Minkowski space-time result in an exponential growth of the amplitudes for long-wavelength fluctuations. 
This growth makes the backreaction important after a certain time eventually shutting-off the instabilities. That is, 
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the non-linearity shut-off the spinodal growth of the modes) 23 |28|. Therefore a non-perturbative treatment of the 
dynamics is required. The large N limit of the scalar sector allows such a systematic treatment, furthermore it is 
renormalizable and maintains the conservation laws ]32]| . 

The gauge contributions to the self-energy of the scalar field can be separated into the contributions from hard 
momentum k sa T gauge modes and those from the soft modes with momenta k -C T '. The hard momentum modes 
will remain in local thermodynamic equilibrium and their contribution to the self-energy of the scalar field to lowest 
order results in a thermal mass eT/\/8[^4j [53|. This contribution adds to that from the hard scalar modes k w T in 
the scalar tadpole diagram that also lead s to a thermal mass squared AT 2 /24. Both contributions are accounted for 
in the thermal mass mf, in equation fl5.l| ) . 

Furthermore, resumming the contribution from the hard modes (in equilibrium) to the longitudinal photon propaga- 
tor (namely the Dyson series) leads to the screening of the instantaneous Coulomb interaction with a Debye screening 
length £d ~ (eT) _1 |Q. This results in the screening of long range forces which would otherwise prevent large charge 
fluctuations. 

The contribution from soft modes of the gauge field to the scalar self-energy will only include non-equilibrium 
corrections through the scalar loops in the photon propagator in the self energy. Therefore the back-reaction of the 
non-equilibrium dynamics of the gauge field onto the dynamics of the scalar field will appear at the two-loop level, 
namely at order a 2 . Therefore we will neglect the non-equilibrium back-reaction of the gauge field upon the dynamics 
of the scalar field. 

The dynamics in leading order in N already reveals the important non-equilibrium features of the evolution. 
Furthermore, we will neglect the back reaction of the vector field on the dynamics of the scalar field since these effects 
are suppressed by at least one power of the fine structure constant a and are subleading in the large N limit. The 
dynamics of the scalar field in leading order in TV is presented in refs. |^3|, 28|, pl[ . 

Since symmetry breaking is chosen along the direction of the neutral field er, we write 



a(x,t) = ^Nip(t)+ X (x,t) ; ( X (x,t))=0 



(5.3) 



where the expectation value is taken in the time evolved density matrix or initial state. The leading order in the large 
N limit is obtained either by introducing an auxiliary field |32| and establishing the saddle point or equivalently by 
the factorizations [E8| 
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The non-linear terms of the a field lead to contributions of 0(1/ N) in the large N limit, and to leading order 
the dynamics is completely determined by the N complex scalars $ r . This factorization that leads to the leading 
contribution in the large N limit makes the Lagrangian for the scalar fields quadratic (in the absence of the gauge 
coupling) at the expense of a self-consistent condition: thus <j> acquire a self-consistent time dependent mass. 
It is convenient to introduce the mode expansion of the charged fields 



$ r (t,2) = 



d 3 k 



a r (k)f k (t) e rk - s + bt(k) f* k (t) e 



(5.4) 



d 3 k 



b r (k)f k (t) e ik - s + at(k)f* k (t)e 



-ik-x 



The equation of motion for the expectation value ip(t) for t > [see eq.(p3.3|)] is given byp3|, 



(5.5) 



(5.6) 



In leading order in the large N limit, the Heisenberg equations of motion for the charged fields translate into the 
following equations of motion for the mode functions for t > |2f| 



fk(t) = 



(5.7) 



The initial conditions for the mode functions are chosen to describe particles of mass M, namely 

f 



A(o) 



A(0) = -iW k f k (0) , W k = y/k 2 + M 2 ; 



(5.8) 
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where M 2 is given in eq.(|5. 

For the vacuum case the initial state is annihilated by a r (k); b r (k) whereas for the initial state of local thermody- 
namic equilibrium at high temperature 



(al(k)a s (k)) = (bl(k)b s (k)) = S rs n k ; n k = 



e t 



With this choice of the initial state we find 



2T 



(5.9) 



(5.10) 



with W k given in eq. ( |5.S| ) and the vacuum case is obtained by setting T = above. 

This expectation value is ultraviolet divergent, it features quadratic and logarithmic divergences in terms of a 
momentum cutoff. The divergences are absorbed in a renormalization of the mass term — > n 2 R and into a renor- 
malization of the scalar coupling A — > Ar. These aspects are not relevant for the discussion here and we refer to 
refs.[^3|, for details. 

After renormalization, ■£? (<j>t$) is subtracted twice, and after rescaling of variables introduced above it is replaced 
in the equations of motion by p3, |28| 



A A f°° f 

— <$t$) ^ AH E( t ) = _A J q 2 dg) | /g(t )|2 coth 



2T 



Q(g-^) 

2q 3 



-^ 2 (0)+^(t) 2 +A fl E(i) 



(5.11) 



and the mass and coupling are replaced by their renormalized counterparts. In order to avoid cluttering of notation 
we now drop the subscript R for renormalized quantities, in what follows /i; A stand for the renormalized qua ntities. 

We now consider the case of a quench from an initial state in which ip(0) = <p{0) = which from eq.(5.6) entails 
tp(t) = as a fixed point of the dynamics J2(], ^8|. This initial case describes either the vacuum case or the case of 
a transition from an initial high temperature phase in local thermodynamic equilibrium. 

After renormalization and in terms of dimensionless quantities, the non-equilibrium dynamics of the charged scalar 
fields is completely determined by the following equations of motion |23, [2q| 



dt 2 
AE(t) 



fJ 2 + q 2 + AE(t) 



/,(*) = ; /,(0) = 



1 



; / 5 (0) = -iW q fg(0) 



q 2 dq{ \f q (t)\ 2 coth 



W a 



2T 



l/,(0)| s 



e(g 2 -M 2 ) 

2g3 



A£(f) 



(5.12) 
(5.13) 



with Wg given by equation (5.2). The full time evolution determined from these equations has been studied in 
refs. ||23 



28 . We summarize the aspects of the non-equilibrium dynamics that are relevant for the generation of 
primordial magnetic fields (for more details see refs.[^2[ §)■ 

In the weak coupling limit A< 1 the back-reaction AS can be neglected as long as E < 1/A. Neglecting the back- 
reaction it becomes clear that the mode functions f q (t) increase exponentially in the band of unstable wavevectors 
q < /i with the result, 



/,(*) 



a„ exp 



(ty/Jj? - q 2 ^ + b q exp (-t^Jn 2 - q 2 ^j 



q<fi 



1 

2JWZ 



1 - 



i W„ 



r 



b q = a* 



(5.14) 
(5.15) 



A feature of the solution (5.15) is that when the exponentially damped solution becomes negligible as compared to 
the exponentially growing one, namely for t > fi , the phases of the mode functions f q (t) freeze, i.e. become constant 
in time and are a slowly varying function of q for long wavelengths, while the amplitude grows exponentially, namely 
the long-wavelength modes behave as 



hit) 



i5 ° e pt 



for q <C (J, 



(5.16) 



This feature of the long-wavelength mode functions will play an important role in the discussion below. 



10 



The exponential growth of modes in the unstable spinodal band q < fi make the back reaction term AX to grow 



and eventually compensate and cancel the 
scale defined by p8| 



-fi in the equations of motion. This will happen at a nonlinear time 



AX(i„; 



(5.17) 



Two important aspects are described by t n i: i) at this time scale the phase transition is almost complete since 
XT,(t n i) = means that \(&$) /2N = fi 2 , namely the mean square root fluctuations in the scalar field probe the 
manifold of minima of the potential. 

ii) At t ~ t n i the amplitude of the field is of order /j/yA and the non-linearities become very important. The back 
reaction AX(i) becomes comparable to and the instabilities shut-off. Thus for < t < t n i the dynamics is described 
by the linear spinodal instabilities while for t > f„; a full non-linear treatment of the evolution is required. However, 
detailed analytical and numerical work both in Minkowski and cosmological space times |^, ^7], have shown that 
the main features of the non-equilibrium dynamics can be gleaned from the intermediate time regime that can be 
studied within the linear approximation. 



Using the approximation (5.16) for the mode functions in the intermediate time regime t n i > t » fi and for 



long wavelengths which give the largest spinodal growth, we find 



AE(t, 



M 2 



32tt 2 M JJi 



2 M t„ 



A for 



'96A for high temperature 



(5.18) 



where we have used n q m T jW q and the long- wavelength approximation W q ~ Wo — y/X/24 T for the high tem- 
perature case. This leads to the following estimate for the nonlinear time to leading logarithmic order in the weak 
coupling A 




32(?r/2) 3/2 

4 V M > ix } 

32(>/2) 3 / 2 
V6A(^ + f ) 



O log log 



(5.19) 



for vacuum 

for high temperature 

The amplitude of the long- wavelength modes at the nonlinear time by the end of the phase transition is approximately 



\f q {tnl)\ - ^= (5.20) 

As we will discuss in detail below this non-perturbative scale will ultimately determine the strength of the magnetic 
fields generated during the phase transition. 

During the intermediate time regime the equal time correlation function is approximately given by 

^\{t)%{t))^\a q \ 2 e 2 ^e-^ . (5.21) 
and its Fourier transform for long wavelengths is of the form 

G > (t,t,x) w |d | 2 e 2/lt (5.22) 
which determines the time dependent correlation length of the scalar field 

((t) - J- . (5.23) 
V A* 

The detailed analysis of the dynamics in refs. |^3| ^8) and the discussion of the main features presented above can be 
summarized as follows: 

• At intermediate times < t < t n \ ~ In the mode functions grow exponentially for modes in the spinodally 
unstable band q < \x. The phase of these mode functions freezes, namely, becomes independent of time and 
slowly varying with momentum. 
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• At a time scale determined by the nonlinear time the back-reaction shuts off the instabilities and the phase 
transition is almost complete. This can be understood from the following argument: the back reaction becomes 
comparable with the tree- level term (for t > 0) when -^{^Q} w /i 2 . This relation determines that the mean 
square root fluctuation of the scalar field probes the minima of the tree level potential. 



During the spinodal stage the correlation length of the scalar field grows in time and is given by equation (5.23). 
This is interpreted as the formation of correlated domains that grow in time, and is the hallmark of the process 
of phase separation and ordering. This correlation length will be important in the analysis of the correlation of 
magnetic fields later. 

The large fluctuations associated with the growth of spinodally unstable modes of the charged fields will lead to 
current fluctuations which in turn will lead to the generation magnetic fields. Thus the most important aspect 
of the non-equilibrium dynamics of the charged fields during the phase transition is that large fluctuations of 
the charged fields associated with the spinodal instabilities will lead to the generation of magnetic fields. Since 
the modes with longer wavelength are the most unstable the magnetic field generated through the process of 
phase separation will be of long wavelength. Furthermore, we expect that the magnetic field generated by these 
non-equilibrium processes will be correlated on length scales of the same order as that of the charged field 
above. 



VI. MAGNETIC FIELD SPECTRUM 



Having discussed the mechanism of generation of magnetic fields through phase separation, the next task is to 
provide a framework to calculate the magnetic field generated. Notice that the expectation value of the magnetic field 
in the non-equilibrium density matrix must vanish by translational and rotational invariance, 

< B\t,x) > p =0 (6.1) 

Hence the relevant quantity to focus on is the equal time correlation function 

C(t, x) =< £*(t, x)B'(t, 0) > p , (6.2) 

where we sum on repeated indices. B{t,x) above is a Heisenberg operator and the expectation value is in the initial 
density matrix. Since the coincidence limit of correlation functions of operators must be defined carefully we define 
the spectrum of the magnetic field as 

S B (t,k) = ~ lim f d 3 x < {B l (t 7 x) 7 B l (t',6)} > p e lts , (6.3) 
z t • >t j 

where { , } denotes the anticommutator. 

In thermal equilibrium the density matrix associated with the time-independent hamiltonian of the system takes 
the form p = e~@ H /Z which is invariant under time translations. Hence, the spectrum S B {t,k) — S B (k) does not 
depend on time and can be computed using the methods of thermal field theory. We will return to this case below to 
provide a consistency check of the formulation out of equilibrium. 

In terms of the transverse component of the gauge field A l T (t,k) — V 13 (k) A J (t,k) with the transverse projection 
tensor V %3 (k) = Sij — the spectrum S B (t, can be written as (with an implicit sum over indices) 



S B (t,k) = ~ I d 3 xe lkS k 2 V%(t,t';x)\.., (6.4) 



2 

in terms of the symmetric correlator of the transverse gauge field 

V%(t,t';x) = V ij >(t,t';x)+V ij <(t,t';x) 
V %3 > (t,t';x) = % < At r (t,x)A 3 T (t',0) > ; V l] < (t,t';x) = i < A 3 T {t! , Q)A l T {t, x) > (6.5) 



Where the correlation functions in (6.5) are Wightmann functions without time ordering computed in the initial 
density matrix. The main reason for introducing the definitions above is that there is a well established framework 
for obtaining these correlation functions in non-equilibrium quantum field theory as discussed below. 
From Ss(t, k) we can extract the magnetic energy density 

/d 3 k 
j^S B {t,k), (6.6) 
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where the ultraviolet behavior is understood to be regulated in some gauge invariant manner, for example with 
dimensional regularization. 

From the magnetic energy density we can define an effective magnetic field B e f f (t) such that 

±Bl f/ (t)=p B (t) (6.7) 

In the non-interacting, thermal equilibrium case, which will be relevant to compare the energy density in the generated 
magnetic field to that in the thermal radiation background, 

S ( °\k) =k [l + 2n(fc)] (6.8) 

where n(k) is the Bose-Einstein distribution function. Subtracting the irrelevant vacuum contribution, the effective 
magnetic field is then given by the Stefan-Boltzmann form 

The phase transition generates a dynamic effective magnetic field AB e ff(t) through the interaction between the 
charged fields and the electromagnetic field. Hence a priori we would expect that AB e ff(t) can be obtained systemati- 
cally in a power series expansion in the electromagnetic coupling a, namely AB e ff(t) — aAB^(f)+a 2 AB^(f) + - • • . 

This would be the case were it not for the fact that the DC conductivity is important to estimate reliably the 
amplitude and the correlations of the generated magnetic field. As it will be discussed below in a high temperature 
plasma the conductivity is determined by fluctuations of the charged fields with typical momenta p ~ T. These short 
wavelength modes are in LTE while only low momentum modes of the charged fields will fall out of equilibrium and 
undergo spinodal instabilities during the phase transition. The DC conductivity is nonperturbative and non-analytic 
in a and to leading logarithmic order it is approximately given bypo], |36| 

CNT 

(6.10) 



a In -^rf 

aN 



with N the number of charged fields and C a numerical constant of order one|35|, |36| . 

Thus while we cannot provide an expression for the generated magnetic field as a power series expansion in a 
because the presence of the conductivity prevents such expansion, the strategy that we pursue here is to treat 
the long-wavelength non-equilibrium current fluctuations in perturbation theory in a while the short wavelength 
contribution will be accounted for in the conductivity. This will be discussed below in detail, but the important point 
of this discussion is that the long-wavelength fluctuations that lead to the generation of the magnetic field will be 
treated to lowest order in a and to leading order in the large N limit. 

The reliability of this expansion will be guaranteed if 

AB ef f(t) 

B eff 

In cosmology the important quantity is the energy density in long-wavelength magnetic fields on scales L equal to or 
larger than the galactic scales. Assuming rotational invariance, we introduce the energy density of the magnetic field 
generated by the non-equilibrium fluctuations 

A PB {L) = — k 2 AS B (k,t) dk . (6.12) 

2tH Jo 

The quantity of cosmological relevance is 

= %(L) 
Pi 

where p 1 = 7r 2 T 4 /15 is the energy density in the thermal equilibrium background of photons. 

As discussed in section IV, the physical magnetic field in a cosmological background is diluted by the expansion as 
B ~ 1/C 2 (n) w ith C (rj) the scale factor. Therefore, in the absence of processes that generate or dissipate magnetic 
fields the ratio fl6.13| ) would be constant under cosmological expansion because of the redshift of the temperature T ~ 
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1/(7(77). In a cosmological background the ratio fl6.13| ) will only depend on time through the generation or dissipative 
mechanisms (such as magnetic diffusion in a conducting plasma) but not through the cosmological expansion. 

For r > 10 -32 the linear (kinetic) dynamo may be sufficient to amplify cosmological seed magnetic fields and for 
r > 10~ 8 the collapse of protogalaxies with constant magnetic flux may be sufficient to amplify the seed magnetic 
fields|, ||, |. 

In the cosmological setting this ratio is approximately constant if the conductivity in the plasma is very large after 
the non-equilibrium dynamics has taken place. The approximations invoked to estimate the spectrum of the generated 
magnetic field will be reliable provided r(L) <C 1. 



A. Non-equilibrium dynamics of electromagnetic fluctuations 

The formulation of quantum field theory out of equilibrium is by now well established in the literature and we 
refer the reader to refs |28|, Q for details. The generating functional of real-time non-equilibrium correlation 
functions can be written as a path integral along a contour in (complex) time. The forward and backward branches 
of this contour represent the time evolution forward and backward as befits the time evolution of an initial density 
matrix p7, |28| , p4|. Gathering all bosonic (scalar and vector) fields generically in a multiplet the generating 
functional of correlation functions is given by [^7[ |34[ 

Z(J+, J_) =y"[P*+P*_] e * J «***{£(*+)-£(*-)+'+*+-'-*-} . (6.14) 

where functional derivatives with respect to the sources J± lead to the non-equilibrium real time correlation functions. 
The doubling of fields with labels ± is a consequence of the fact that the non-equilibrium generating functional 
corresponds to forward (+) and backward (— ) time evolution and suggests introducing a compact notation in terms 
of the doublet and the metric in internal space Q 

c ab = (j \) = c-J , (6.15) 

where the labels a, b = ±. This notation is particularly useful to obtain the non-equilibrium version of the Schwinger- 
Dyson equations for the propagators . 

In the case under consideration the main ingredients for the non-equilibrium description are the following 

• Transverse photon propagators 

The real time Green's functions for transverse photons are given by 

(4^(z,*)Ag0rV)> = -if (0 e-*(*-**'> , (6.16) 

where the explicit form of T>fj(k;t,tf) is 

V±+(k;t,t')=V ij (k) [2?*(M,Oe(*-0+^(M,«W-*)] > ( 6 - 17 ) 
VT j -(k;t,t')=V ij (k) [D>(k;t,t')e(t' -t) + V<(k;t,t , )e(t-t')] , (6.18) 
V+-(k;t,t') - Pijtf) V<(k;t,t') ; V7. + (k;t,t') = V l3 (k) V>{k;t,t') (6.19) 

and Vij (k) is the transverse projection operator, 

V ij (k) = S ij -!^.. (6.20) 



• Scalar propagators: 



1 H^k - 

^(x,t)^ b Hx',l/)) = -tS rs I ^3 G£{t,t') e -*(*-"0 , (6.21) 
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= g> (t, t')e(< - <') + g< (t, - t) , 


(6.22) 


G- k -{t,t') 


= G>(t,Oe(t'-t) + G<(M')e(t-0 - 


(6.23) 


Gfc"(*,* / ) 


= G<(t,t') ; G-+(t,t') = G>(t,t'), 


(6.24) 


(*,f) = 


^{[l + nd/fcCO^CO+nk/^O/fcC*)} , 


(6.25) 


G<(t,t') = 


= 5 {[1 + n fc ] + n h f k (t)fffl)} • 


(6.26) 



For the scalar propagators we h ave used the expansion in terms of mode functions given by eq. (5.5) and the 
expectation values given by eq. (5.9). 



The scalar propagators given above imply a non-perturbative sum of cactus-type diagrams when the expectation value 
of the a field vanishes. These propagators are depicted in fig. [|. 




FIG. 1: Scalar propagators in leading order in the large N. For vanishing expectation value of a these are given by a sum of 
cactus-type diagrams. 



B. Non-equilibrium Schwinger-Dyson equations 

In this section we derive an exact expression for the spectrum of the magnetic field from the non-equilibrium set 
of Schwinger-Dyson equations. 

The full non-equilibrium propagator for the photon field is obtained from the non-equilibrium effective action 
resulting from integrating out the charged fields, and which up to quadratic order in the photon field is given by 



TL4 d 



d^xd^y {Al T {x) [d^5 A {x 



y)c a % 3 



Il% T (x,y)]A b jtT 



(y)} 



(6.27) 



with an implicit sum over all repeated indices. 

To lowest order in a and to leading order in the large N expansion, the non-equilibrium contribution from the 
scalar fields to the photon polarization is given by the two diagrams shown in fig. |^. 

The tadpole diagram (a) in fig.0 gives a contribution that is local in time, thus we write the polarization 



n 



/ ad 



(t) c ab 8fj 5\x 



•y) + IE? T (a? 



y°) 



(6.28) 



where the form of the local (tadpole) contribution H tad (t) in terms of the metric c ab is necessary for consistency 
with the form of the effective action since the kinetic term is proportional to the metric |j4j and we have used spatial 
translational invariance. The contribution II is non-local, and in equilibrium it fea tures a bsorp tive p arts. T he di fferent 
components a, b = ± are defined in the same manner as the Green's functions J6l7|)-(|sl9j) and (^2§)-(|jg). The 
contribution to order a to the non-loca l par t of the polarization is depicted in diagram (b) in fig.§. 

The photon propagator given in eq. (6.16) is the inverse of the operator in the quadratic effective action above and 
obeys the non-equilibrium version of the Schwinger-Dyson equation [j34| 



d 4 y [d^S 4 (x~ y) c ab d i: 



V$(x,y)] Vl k c (y,z) 



6 4 (x - z) 



(6.29) 
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(b) 



FIG. 2: Photon polarization to lowest order in a and to leading order in the large N. The scalar loops are in terms of the 
scalar propagators displayed in figjl] 



From the expressions (6.17)-(6.19) for the different components of the photon propagator it proves convenient to 
introduce the following combinations 



V^(x,y)=V>(x,y)-V<(y,x) 
Vf j {x,y)=Vf j {x,y)6{x Q -y Q ) 
Vf j (x,y) = V? j (x,y)6(y°-x ) 
V? j (x,y)=V>(x 1 y)+V<(y,x) 



antisymmetric propagator 
retarded propagator 
advanced propagator 

symmetric propagator. 



(6.30) 
(6.31) 
(6.32) 
(6.33) 



and similar definitions for the polarization II. 

Using spatial translational invariance we can simplify the above form of the Schwinger-Dyson equations by intro- 
ducing the spatial Fourier transform of the transverse propagators 



V% T (t,t';x-x') 



d 3 k 
(2^)3 



jk-{x-x') v V ab (t,t';k) 



and introducing the spatial Fourier transform of the transverse polarization tensors defined by 



Il a T b (t 1 t';x-x , ) = ± 



d 3 k 



(6.34) 



(6.35) 



We can now write down the Schwinger-Dyson equations for the tr ansve rse photon propagators. Since there are only 
two independent functions D >,K the Schwinger-Dyson equations (6.29) can be simplified to a set of two equations 
for a particular combinations of these. Some straightforward algebra leads to the following set of non-equilibrium 
Schwinger-Dyson equations (we have suppressed the subscript T but the equations below refer to the transverse 
components) 



dt\ 



— + k' + n tad {t)) S(t-tx)+U R (t,tv,k) 



V R {tx,t';k)=8{t-t') 



and 



dti 



dtj 



+ k 2 + U tad (t) S(t - h) + U R (t, h;k) 



V H {t u t';k)= / dt x n H (t,t i; fc)X)A(ti,t';fc) 



(6.36) 



(6.37) 
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with the definition 



U a (t,t 1 ;k)= n > (i,ti;fc)-n<(t,t 1 ;fc) 0(ti - i) 



(6.38) 



A remarkable aspect of this set of equations is that the right hand side of eq. (|6.37 ), namely the inhomogeneity in 
the set of equations, only involves the non-local contribution to the polarization, which to lowest order in a is given 
by diagram (b) in fig. |. This is a result of the form of the local (tadpole) contribution which is proportional to the 
metric c ab . This is an important point to which we will come back later in the discussion. 

The antisymmetric propagator T>c(t,t";k) is the odd solution of the homogeneous evolution equation 



dt 2 



n tad (t)^j V c (t, t"; k) + J dt' IL R (t, t', k) V 



c(t',t",k) = 



with the constraint 



dV c (t,f,k) 



01 



1 . 



(6.39) 



(6.40) 



This relation guarantees the correct equal time canonical commutators. 

Since the kernel of the integra l equ ation on the left han d sid e of eq. ( |6.37 ) is the same as for the equation that 
defines the retarded propagator ( 5.36j ), the solution to eq. (6.37) is given by 



V H (t,t';k) = / dh / dt 2 V R (t,h;k)n H (h,t 2 ;k)V A (t 2 ,t';k) + F(t,t';k) 



(6.41) 



where T> B (t,t';k) is the solution of eq. (3.36) and the function F(t,i/;k) symmetric in t and t' is a general solution 
of the homogeneous equation 



dh 



— + k 2 +U tad (t)) 5{t-h)+^ R (t,h;k) 



F(h,t',k) = 



(6.42) 



The homogeneous solution F can be constructed systematically and its physical significance will be discussed below. 

N ow are now in position to provide the final expression for the spectrum of the magnetic field. From the expression 
(5.4) and the solution found above, we find 



S B (t,k) = -ik 2 V H {t,t';k)\ 



(6.43) 



The expression ( |6.41 ) can be simplified further by taking into account the theta functions in the definitions of the 
retarded and the advanced (T>r and T>a) propagators [see eqs. ( |6 .31 )-( |5~3"2 )] as well as the antisymmetry T>c(t, t'; k) — 
—T>c{t', t; k) of T>c leading to 



S B (t,k)=i k 2 



( dh j dt 2 Vc(t,t2;k)Vc{t,h;k)U H (h,t2;k)- F(t,t;k)\ 

J to J to ) 



(6.44) 



Since the product T>c{t 1 t2]k)T>c{t,h',k) is symmetric in the exchange h h we can replace n#(ti,t2; k) by 
2n>(ti, ti\ k) and write the final form of the spectrum separating, for further convenience, the contribution from the 
inhomogeneous and homogeneous solutions to Sb- 

S B (t,k) = S I B (t,k) + S B r {t,k) 

S B (t,k) = 2ik 2 / dh dt 2 Vc{t,t 2 ;k)Vc{t,h;k)Tl > (h,t2;k) 

J to J to 

5f (t, k) = -ik 2 F(t, t; k) . (6.45) 

where £o is some initial time before the phase transition and n> does not include the local tadpole contributions, it 
is the non-local part of the polarization. 

There is an important aspect associated with the homogeneous solution F(t,t',k) and its contribution to the 
spectrum of generated magnetic fields S^(t, k). This aspect is revealed by noticing that the expectation value of the 
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transverse gauge field Ar(t, k) — {Ax{t, k)), namely the mean field obeys the same homogeneous equation of motion 
as F(t,i/,k), 



dt\ 



2 +k 2 +U tad (t) ) 5(t-tx)+n R (t,tx;k) 



A T (t,k) = 



(6.46) 



(we have suppressed the vector indices to avoid cluttering of notation) . 

Thus the homogen eous solution F(t,t',k) can be constructed out of the independent solutions of the mean field 
equations of motion (6.46). The main reason that we bring up this point is to highlight that the solutions to the 
mean-field equations of motion are only part of the contributions to the generation of magnetic fields through non- 
equilibrium processes. However, as it will be discussed in detail in the next section, this contribution can be neglected 
in the present case of non-equilibrium spinodal decomposition in many circumstances, and the term Sg(t, k) dominates 
for late times. 



C. Electric fields: 



For completeness we now address the generation of electric fields. The importance of generation of electric fields is 
mainly related to the question of equipartition. It is often assumed that the energy density stored in electromagnetic 
fields is equally partitioned between electric and magnetic field components, namely between temporal and spatial 
gradients. While this is usually the situation in equilibrium, it is not necessarily the case strongly out of equilibrium 
such as the situations envisaged in this article. 

The electric field is the Hamiltonian conjugate field to the vector potential, and its transverse component is given 
by Ej,(t, k) = —Aj,(t, k), therefore the equal time correlation function of the electric field is given by 



s ET {t,x) =< Ai r (t,x)Ai r (t,d> c =d t d t/ < >< 



(6.47) 



where just as in the case of the magnetic field we wrote the equal time correlator as the symmetrized connected 
two-point correlation function. 

Following the steps leading to eq. J6.44 ) for the magnetic field, and using the following identities 



d t V R (tM\k) = d t Vc{tM)k) 6{t-t x ), d t ,V A (t 2 ,t';k) = -d t ,V c {t' \t 2 ;k) 9(t 2 - t') 

we now find 

S ET (t,k) =dtdv \li f dh f dt 2 D c (t,t2;k)Dc(t',t 1 ;k)n > {t 1 ,t 2 ;k)- F(t,t',k)X 

I J to J to J t=f 

The number of photons produced through the non-equilibrium processes is given by 



N(t, k) 



S ET {t,k) + S B (t, k) 
2k 



(6.48) 



(6.49) 



after summing over the two polarization states. 

The final form for th e spe ctrum of magnetic and electric fields given by eq. ( 3.45 )-( 5.48] ) as well as the number of 
photons given by eq. ( 3.49 ) are the main tool to compute the spectrum of electromagnetic fields generated during 
non-equilibrium processes and one of the main results of this article. 

We emphasize that these expressions are exact and general and apply (with rather minor modifications as explained 
in a previous section) in the cosmological setting. They allow to study the problem of the generation of magnetic 
fields through non-equilibrium processes in general situations. 

Furthermore, the final expressions for the spectrum of electromagnetic fluctuations ( 6.45| . |6.48 ) is valid more gener- 
ally in spinor electrodynamics since it only involves the full polarizations and the Schwinger-Dyson equations for the 
correlation functions which ultimately lead to the final expressions are general. 



D. Spectrum of fluctuations in equilibrium 



Before focusing on the study of the magnetic field generated during non-equilibrium phase transitions, it is both 
illuminating as well as important as a consistency check to address the case of thermal equilibrium. In this case 
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the simplest manner to compute the spectrum is by using the imaginary time or Matsubara formulation, where the 
transverse photon propagator is written as pij 

f ,73 J, _ t 

V%(t,x) =< A^t,x)A' t (t',6) >= / ——-7 e lk - s T V^ffc) V T (u n ,k) e -*"-(r-r ) ; u n = 2nT n (6.50) 



with 



PrK,*0= / ^0 g (fc '. fco) (6.51) 

and pT{k,ko) is the spectral density for transverse photons which is an odd function of kg. The spectrum for the 
magnetic field is obtained from the equal time limit of the Euclidean propagator, and is therefore given by 

S B (k) = 2 k 2 Tj2v T (uj n ,k) . (6.52) 

nGZ 

The sum over the Matsubara frequencies can be performed using the methods described in |m[ leading to 

S B {k) = 2 [ dui k 2 n(u) pr(u), k) . (6.53) 



The general form of the spectral density in terms of the transverse polarization is given by [ p4| 

/, , s 1 Imn T (fc,fc ) 

PT{k,k ) = 5 - (6.54) 

I" [k 2 -k 2 - ReII T (fc, fc )] + [Imn T (fc, k )} 2 

The spatial and temporal Fourier transform of the retarded and advanced transverse photon propagators are given 
by g 

Vf j (k,uj)=V ij (k)V R {k,uj) ; V R {k,uj) ' 



-co 2 + k 2 + Ren T (fc,w) + i ImIT T (fc,cj) 

V A (k,uj)=V ij (k)V A (k,oj) ; V A (k,oj) = 2 ,, 2 , RTT ,) : t T\ (b ^ ( 6 ' 55 ) 
J — uj a + k l + RclLr(fc, uj) — i imllT(/e, uj) 

We are now in conditions to establish contact with the non-equilibrium result of the first section. In equilibrium the 
polarization and the propagators arc functions of the time differences and we can then take the Fourier transform in 
time of T>H(t, t' , k) given by eq. ( 6.41 ) 

T> H (k,w) = T> R (k,u) [U > {k,uj)+U < {k,uj)]V A {k,uj)+F{k,uj) (6.56) 

where we have explicitly written II in terms of n< >. In equilibrium the detailed balance (or KMS) condition relates 
these components of the polarization to the imaginary part p3j] 

n>(fc, uj) = - [1 + n(uj)} ImII(fc, uj) ; Tl < (k,uj) = - n(u) kall(k, uj) (6.57) 

7T 7T 

with n(uj) being the Bose-Einstein distribution function. 

Furthermore, the Fourier transform of the homogeneous solution F(t, t' , k) obeys the Fourier transform of eq. ( 6.42 ) 
namely 

V^(u,k) F(uj,k) =0 (6.58) 

Combining all the a bove ingredients and using the fact that ImII(fc, uj) is an odd function of we are led to the 

conclusion that eq. ( 6.43| ) for the spectrum becomes 

S B (k) = I dw {2 k 2 n{uj) p T {uj, k) + k 2 F{uj, k)} (6.59) 



with pr(uJ, k) given by eq. (6.54). 

This result differs from that obtained via the equilibrium propagator eq. (6.53) by the contribution of the homoge- 
neous solution F. 
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In order to understand the source of the discrepancy between the two formulat ions, namely the contribution of the 
homogeneous solution F(k,u>), let us focus on the denning equation for F(k,oj) ( |6.58 ). 



In order for a non-vanishing solution to this equation to exist, from the expression for 25r(/c,u;) we infer that 

lo 2 ~k 2 ~ReH T (k,Lu) = (6.60) 
lmH T {k,uj) = (6.61) 

Eq. ( 6.60| ) determines the dispersion relation of quasiparticles and eq. (3.61) determines that these quasiparticles 



must have zero width, i.e, the solut ions of the homogeneous equations are propagating quasiparticles with zero width 
and a dispersion relation given by ( |6~60l) . The homo geneous solution is therefore 

F(k,u) oc 5{uo 2 - k 2 -ReII T (fc,cj)) (6.62) 

In a non-perturbative resummation of the Dyson series for the propagators and the spectral density, the limit when 
Imlly — * 0* leads to 

p T -> tS(lo 2 - k 2 - Ren T (fc, u)) (6.63) 

and we reco gnize in this case that the possible contributions of the homogeneous solutions are accounted for in the 
expressions and arise from the propagating states, and to zeroth order account for the contribution to the 



spectrum of the magnetic field from free photons. 

However, in a perturbative expansion perturbation theory is only reliable away from the single particle poles in the 
propagator and the limit ImllT — > ± will miss the contribution from the isolated poles. The homogeneous solution 
F(k,uj) gives the required contribution, thus guaranteeing the consistency of the perturbative expansion. Thus while 
the homogeneous solution F must be included in a perturbative calculation of the spectrum, it will be accounted for 
in a non-perturbative computation that includes a Schwinger-Dyson resummation for the full propagator and should 
not be included in the computation of the spectrum. 

In a plasma the true degrees of freedom are quasiparticles, in particular for the power spectrum of magnetic 
and transverse electric fields, the important degrees of freedom are transverse plasmons. In the hard thermal loop 
approximation f2~if the transverse plasmons contribute to the spectral density prjk, to) a pole term of the form 
Zt{K)[8{uj — uip(k)) — S(uj + uip(k)} with u> p (k) being the plasmon dispersion relation E4fl with co p (0) = eT/3. For 
k 3> eT, Lo p (k) ~ k ; Zt ~ 1 /2k and we recover the power spectrum of free field theory. However for k < eT the power 
spectrum of the magnetic field reveals the presence of collective degrees of freedom, 5e(fc) ~ kZx(k)[l + 2n(cu p (k)] 
for e 2 T <k<eT. 

For ultrasoft frequency and momenta k, to <C eT which is the case of relevance for the cosmological case, the effective 
low energy form of the gauge field (transverse) propagator (well below the plasmon pole) is 

V(k,u)~^— I — (6.64) 

lo 1 — k z — icruo 

leading to the small frequency, long-wavelength form of the spectral density 

p T (k,U})~-— r^W—f so- ( 6 - 65 ) 

Therefore, for to -C cr 

k 2 n(Lo) PT {k^)^-^-—2 ; r= — . (6.66) 

Thus for long wavelengths k <§; a w T/a the spectral density features a pole at zero frequency leading to a long- 
wavelength power spectrum 

S B (k<^o-)=T, (6.67) 

coming from the first term in eq.( 6.59| ). 



The plasmon pole contributes to S*s(fc) through the homogeneous solution F in eq.(6.59) yielding, 

— - 1 + 2n(w p (0) ~ ~ - 

w p (0) w^(0) e 2 T 



SBplasmon - —JT^i 1 + 2n(w„(0)] - - (6.68) 



where we have used the long-waveleng th lim it of the dispersion relation Lo p (k) and residue (fc) pi)] ■ The contribution 



( 5.68 ) is clearly much smaller than ( 6.67 ) in the long- wavelength limit k -C eT. Hence, the magnetic field power 
spectrum in equilibrium is Ssik) — T. 

Having clarified the role of the homogeneous solution and the form of the power spectrum within the more familiar 
equilibrium setting, we return to the non-equilibrium situation. 
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VII. SPECTRUM OUT OF EQUILIBRIUM 



The set of equations for the spectra of magnetic and electric fields ( |6.45 )-( 5.48 ) is exact and only involves the 
transverse polarization from which the full propagator and inhomogeneous and homogeneous solutions are obtained. 
Obviously in order to make progress and obtain an estimate for magnetic and electric field generation one has to make 
approximations. In what follows we will obtain the spectra for magnetic and electric fields to lowest order in a and 
to leading order in large N. 

To this order the tadpole (local) and bubble (non-local) contribution to the polarization are given respectively by 
diagrams (a) and (b) in fig. |2j. Their explicit expressions are given by 



n tad (t) = -ie 2 N 

II«,>(ti,*a,£) = ie 2 N 
which leads to the transverse components 
n>(fi,t2,fc) = 2ie 2 N 



d 3 q 
(2tt) 3 
d 3 q 



G>(t,t,q) 



(2tt) 3 



(2 qi + k t ) (2 qj + kj) G>(ti,t 2 ,g) G > (h,h,\q + k\) 



(7.1) 



d 3 c 



(2tt) 3 

n<(ti,ta,fe) = n>(t 2 ,ti,fc) ■ 



q 2 (l-cos 2 0) G > (t 1 ,t 2 ,(7) G>(ii,i 2 ,|g + fc|) 



(7.2) 



where cos 9 = q ■ k. The scalar pro pagat or G > (t,t / ) in terms of the mode functions f q (t) that satisfy the evolution 
equations (5/7) is given by [see eq. ( 6.25 )], 



G>(ii,i 2 ;fc) - - [(l + n fc )A(*i)/ fe *(* 2 ) +n k ft{tx)h{t 2 )] 



(7.3) 



and we assumed that before the phase transition the scalar fields have occupation rife (obviously this assumption can 
be relaxed straightforwardly). 

From these expressions we can obtain a reliable estimate of the non-equilibrium effects in the polarization. For 
intermediate times after the phase transition, when the dynamics is dominated by the spinodal instabilities and before 
the non-line aritie s in the scalar field evolution are important, the long- wavelength mode functions are approximately 
given by eq.( 5.16| ). 

Near the end of the phase transition for t ~ t n \ w In ( jj the leading order time dependence of the scalar Green's 
functions is then approximately given by 



G>M',fe)t,t'~t„ 



A 



(7.4) 



which allows us to estimate the order of magnitude of the different terms in the polarization. For the tadpole 
contribution (local term) we find 



W ad {t) t 



e 2 (i 2 e 2 ^ ^e 2 ^- + 0( — 



1 

7x 



(7.5) 



This estimate is consistent with the fact that the tadpole contribution is e 2 < > and near the end of the phase 
transition the mean square root fluctuations of the scalar field probe the vacuum state, namely < $ 2 >^ /i 2 /A. 
Similarly, we find for k -C /J, 



n> (t, t' , k)t~t„i — 



128 7T2 



M 2 
M 2 " 



;(t+t')(2/ 1 -L) 



l + O 



(7.6) 



which for t ~ t' ~ t n i is of the order 1/A 2 . 

Since t he re tarded polarization given by eq. ( 6.38j ) is antisymmetric in the time arguments the leading contribution 
given by (7.6) above, actually cancels since II<(t, t',k) = H > (t',t, k). A detailed analysis (see ref. jlj|) reveals that 
the retarded non-local part II/j is of the same order as the tadpole contribution, namely 



.2 t_ 

' A 



(7.7) 



and that in fact at long times there is an exact cancellation between the two terms for k « in Minkowski space-time. 
This is in agreement with the equilibrium result of a vanishing magnetic mass urn . 
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A. Spinodal decomposition from a vacuum state 



To understand the order of magnitude of the generated magnetic fields and to gain insight into the main features 
of the non-equilibrium processes, we now study the situation of spinodal decomposition from an initial vacuum state. 
This corresponds to considering an initial state in which the a field rolls down the potential hill while the fluctuations 
of the charged fields are in the vacuum state at the initial time. The dynamics of the scalar field in this case has been 
studied thoroughly analytically and numerically in references [ p0[ p3[ p7[ p8| to which the reader is referred. 

F or ti mes s maller than or of the order of the nonlinear time (5.19) we can evaluate the tadpole diagram (7.1) using 
eqs. Q and fjj) with the result, 



W ad (t) 



e 2 N 
64 7ri 



(p 2 +M 2 ) 



l + 2n 3 

ti (Aft) I 



1 + 



Let us now consider the retarded bubble diagram 

d 3 q 



n fl (ti,t 2 ;fc) 



-Ae 2 N 



(2tt) 3 



q (l — cos 



111! 



G>(ti,t2,g) G > (t 1 ,t 2 ,\q + k\) 0(t 2 -h) 



(7- 



(7.9) 



where we used eq.( p.3S| ) and (|7.2[). 

We are interested in the generation of long-wavelength magnetic fields with k <C (m. The bubble diagram n^(ti , t 2 ; k) 
takes the form, 



n fl (ti,*2;*) = -4 



+oo 4 



q dq 



+i 



2tt) 2 
•° qUq '+ 1 



dx{l-x 2 )lm G > (t 1 ,t 2 ,q)G > (t 1 ,t 2 ,\q + k\) 0(t 2 - ti) 



l 7^ J_ x dX (1 " } Im i fq{tl) f ^ f W+^ f W+Z\M 9 ^ 2 - *0 • ( 7 - 10 ) 

where x = q ■ k and we set the temperature equal to zero. For times smaller or of the order of the nonlinear time 
(5.19) we can evaluate this bubble diagram using eq.(5.15) for the mode functions. This gives for pit ^> 1, fixed t' 
and k <C fJ., 



Im 



fS) hV)f w+ 0) J w+ %\(t') 



= e x ^ ' 



2 \a 2 q \ lm(a 2 ) e^fe+T*) + Im(a*) e 



Calculating now the g integrals in eq. ( |7T0| ) for fi t > 1 yields 
n fl (t,t';fe)'" = " ■<- V 



16 Af ^/7j7T2 ts 



1 + 2e -t'( 2il -fc) 



(7.11) 



where we use d eq .( ^Tl ). 



Using eqs. (7. 8) and (7.11), the photon evolution equation (3.39) takes the form 



j2 2 tit 

^ + k 2 + e 2 C 26 -^ 
at 2 t2 



V c (t,t";k) 



*(2M-fc) 



t2 



e 2 L(t,t") = , 



(7.12) 



where 



C 2 = 



N n 2 +M 2 



L(t,t") 



64 7T^ M y/ji 

N(n 2 + M 2 ) 



t" 



dt' 



1 + 2e 



-f(2/*-fc) 



V c (t',t";k) 



(7.13) 



Eq.(7.12) is valid for times larger than l//x and smaller or approaching the nonlinear time t n i [see eq.(5.19)]. It can be 
considered an inhomogeneous equation for T>c{t' , t"; k), where the inhomogeneity is to be determined self-consistently. 
Let us first consider the homogeneous version of equation ( [7.12 ) , 



d 2 )2 
d^ + k 



e 2 C 2 



=2 tit 

IT 



4>{t) = o 



(7.14) 
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This equation has as general solution for // t 3> 1 and k <C /i 

0(0" 



3 H * 










+Bsta CS e ")] 









(7.15) 



where A and i? are ar bitrar y constants. 

In order to solve eq.( 7.12| ) we use the Green function of eq.(7.14) 



2jxt 



4 + * 2 + e 2 C 2 ^ 
dt 1 t2 



G(t, t') = S(t - t') 



that can be written explicitly as, 



6 O 



sm 



,m4 



COS 



eC 



(7.16) 



We can thus write the antisymmetric solution of eq.(7.12) as follows 
1 



V c {t,t"-k) = (tt')l e~f ('+*') sin 
e C 



eC ( e^ 1 e^ 1 *" 

"7*" vTT ~ 7"F 



eft' 



G(M') 



t'(2M-fc) 



L(t', t") - (t «-> t") 



_ (7-17) 

The first term provides an antisymmetric s olutio n of eq.(7.14) which is normalized according to eq.(3.4C). The second 
term comes from the inhomogeneity in eq.(7.12). 

For t ~ i n ; the first term is of the order A*. Let us now compute the second term which w ill do minate in such 
regime. From the explicit form of L(t' , t") ( [7.13 ) and the solution for T>c(t, t"; k) given by eqn. ( 7.17 ) it is clear that 
L(t',t") is a slowly varying function of its arguments. In particular we can set t' = t n i and take L(t',t") ~ L(t n i,t") 
outside of the integral. 

Using the results given in eqs.(7.16) and (7.17) we find the following result, 

^2 at' 



dt' G(t,t') 



e *f 



c 2 t + ^ 



(7.18) 



for l//i <C t < t n i . We find at t = t n i , 



dt' G{t nU t') — 



2p,t' 



e 4 ^ 



2e 2 C 2 t nl ^ 



(7.19) 



where we have neglected rapidly oscillating phases. 

We finally get collecting all factors and taking into account that L(t' ', t") varies slowly with t' compared with e 2 ' 1 * , 



T> c {t,t ;k) = tt^- e 



4n 



C 2 t 1+ ^ 



and 



Vc(t nl ,t";k)= L{tnl f 3 \ 2 
2 C 2 i nI ^ 



(7.20) 



(7.21) 



Notice that e 2 cance ls in eqs.( 7.20 ) and ( 7.21 ) since both integrals ( 7.18| ) and ( 7.19 ) are proportional to e -2 and the 
inhomogeneity in eq.(7.12) is proportional to e 2 . Therefore, T>fj(t,t', k) is of the order e° for 1/fi <C t < t n i, k <C n 
and fixed t' . 

We thus obtain for the magnetic field spectrum in the long wavelength limit from eq.( |6.45 ), 



S B (tm,k) k ^ 2ik 2 



dti 



dt 2 v c (t n i ,h;k)T> c (t n i , h ; k) n> (ti , t 2 ; k) 



(7.22) 
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and w e use d the fact that the homo gene ous contribution Sg (t, k) ineq.(^45) can be neglected for t = t n i accor ding 
to eq.( 7.15 ). We can now insert eq.( 7.21 ) for T>c(t n i, t"; k) into eq.( 7.22| ). The fastly growing factors e 2tLt in eq.(7.6) 
combine with powers of t to give near t n i ~ ^ log ^ a factor 



A 2 (logi) 2 + ^ 

dominating the integrals. We used here that L(t n i, t) ~ (t n i — t) 2 for t ~ t n i. This relation follows from the fact that 

T>c(t, f; k) t =*' (t ~ f) due to eq (|6.40| ) and using the relation (|7.13| ) for < -> t". 

Performing the integrals in cq. ( 7.22[ ) for t = t n i we obtain neglecting corrections in 1/ In j , 



SB(tnl,k) 



Ne 2 



X 2 (logi) 2+ ^ 



2 9 7T ' 



1 + 



/ 1 



Vlogi 



(7.23) 



where we used eq.( 7.13 ) . 

We read in eq.(7.23) the correlation length of the scalar field at the nonlinear timepC 




(7.24) 



The expression ( 7.23Q for the spectrum of the magnetic field clearly reveals that the correlation length of the magnetic 
field tracks that of the scalar field during the stage of spinodal decomposition. This result is intuitively clear since the 
magnetic field is generated by the long-wavelength spinodal instabilities corresponding to the formation of correlated 
domains that grow in time during the spinodal regime as ~ \j2tj [i P0(| . From Ssjt , k) we can extract the 
magnetic energy density in long- wavelength scales equa l to or larger than L through eq.( 3.12] ). 

The energy density on scales > L given by eq. ( 6.1 2| ) can be computed in closed form in the limits L ^> £(t n ;) or 
1/fx < L < \{t rd ) using eq.(^2|). We find 



Pb{L) 
Pi 



Ne 2 



X 2 log 4 I 2 11 y/2 5tt 5 



and 



Pb{L) 
Pi 



Ne 2 



X 2 log ^ 



I 32 

A 



M 2 



PL (LT) 



(7.25) 



(7.26) 



for the energy on macroscopic scales . This latter can be interpreted as energy on magnetic fields coherent on scales 
of the order or bigger than L. We see that large scale magnetic fields are strongly suppressed as L~ 5 . 



We now compute the spectrum of the electric field from eq.( |6.48 ). We evaluate the integrals is eq.(6.48) by the 
same lines as the eq.(7.22) using eqs.( 7.2fj| ) and ( [7^ ) for Dc(t,if;k) and H>(ti,t2] k), respectively. We find that the 
end-point near f„; dominates the integrals with the result 



Se t (0,tni) 



N e 2 n 



log! 



(7.27) 



and therefore, 



S B {k,t nl ) _ k 2 1 
SE T {k,t nt ) [i 2 log A 



(7.28) 



Thus we see that for long-wavelengths k <C [i the strength of the electric fields generated are much larger than those of 
the magnetic fields. This result clearly indicates a violation of equipartition as a consequence of the non-equilibrium 
generation of electromagnetic fields. If spinodal decomposition occurs from a vacuum state, namely the a field rolls 
from the top of the potential and the charged fields are in a vacuum state, the non-equilibrium processes generate 
mainly electric photons. 
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B. Spinodal decomposition in a high temperature plasma 



We now study the case of a rapid (quenched) phase transition from an initial high temperature phase to a final low 
temperature phase. Before the phase transition the system is in (local) thermodynamic equilibrium at a temperature 
T » T c with T c oc \xj\p2A\ and the long- wavelength k <^T modes acquire a thermal mass 0] 

m T = \[^T. (7.29) 
The occupation number of a mode of momentum q is given by the Bose-Einstein distribution function 

n(q) = ; w(g) = Jq 2 +m 2 T (7.30) 

e t — 1 

Short wavelength modes with momenta q » XT remain in local thermodynamic equilibrium during the (quenched) 
phase transition, while long- wavelength modes undergo critical slowing down pl[ , freeze out and fall out of equilibrium. 
These long- wavelength modes will undergo spinodal decomposition and their mode functions for pT 1 < t < t n i will 



grow as (see eq. (5.15) 



f q (t)^-ia q e^-<?/^\ aq ^ a ^^±. (7.31) 

since for T 2 > T 2 = ^[i 2 it follows that m T > fx. 

In a high temperature plasma, charge fluctuations lead to a large conductivity in the medium. In equilibrium 
the conductivity is obtained from the imaginary part of the photon polarization and it is dominated by particles of 
momenta p ~ T in the loop with exchange of photons of momenta eT < k <C T |35|, [36) . The Drude conductivity for an 
ultrarelativistic plasma at temperature T much larger than the mass of the charged fields (in this case itlt ~ 
is given by 

e 2 N(T) t 

(7.32) 



with J\f(T) oc T 3 is the number of particles plus antiparticles in the plasma and r is the transport relaxation time. A 
naive estimate based on the scattering amplitude with the exchange of photons or charged particles with momenta 
p ~ T would indicate that t _1 ~ JV(T) a 2 /T 2 ~ a 2 T (with a 2 /T 2 the typical cross section from the exchange of a 
particle with p ~ T) and would lead to a Drude conductivity of the form 

a ~ - (7.33) 
a 

However, a careful analysis including Debye (electric) and dynamical (magnetic) screening via Landau damping leads 
to the conclusion that the conductivity is given by J3^, |3(]] 

C NT , 

°=^m (734> 

with N the number of charged fields and C ~ 0(1). Thus we consider separately the contributions to the photon 
polarization from loop momenta in the two very different regimes: a) the hard momenta p ~ T 3> /i correspond 
to charge fluctuations that are always in local thermodynamic equilibrium, b) the soft momenta p <C /i fall out of 
equilibrium and undergo long- wavelength spinodal instabilities. The contribution from hard momenta will lead to 
a large equilibrium conductivity in the medium, while the contribution to the polarization from soft momenta will 
contain all the non-equilibrium dynamics that lead to the generation of electromagnetic field fluctuations. 

As the instabilities during the phase transition develop, the fluctuations of the charged fields will generate non- 
equilibrium fluctuations in the long-wavelength components of the electric and magnetic fields and the ensuing gen- 
eration of long-wavelength magnetic fields. However, the large conductivity of the medium will hinder the generation 
of electromagnetic fluctuations, hence the conductivity must be fully taken into account to assess the spectrum of the 
magnetic and electric fields generated during the non-equilibrium stage. 

In equilibrium the long- wavelength and low frequency limit (k, uj — * 0) of the spatial and temporal Fourier transform 
of the transverse polarization is given by 



Tlrik, w) — iujo 



(7.35) 
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Thus we write for the full transverse polarization for long-wavelength electromagnetic fields 



n r (t, t', k) = a —5{t - 1') + n noneq (t, t', k) 



(7.36) 



with TLnoneq(t, f , k) the contribution from the spinodally unstable long- wavelength modes given by equations ( |7.l| ). 
Our strategy is to obtain the non-equilibrium contribution to the spectrum of electromagnetic fields to lowest order 
in a but treating the conductivity exactly. 

The zeroth order propagators are now obtained by considering the contribution from the conductivity to the 
equations of motion, the relevant retarded, advanced and symmetric transverse correlators obey 



r d 2 (2 d 
r d 2 l2 d 

d£ +k+a dt 



+ fc + a— 



V ( ° ) (t,t',k) = S(t-t') ; V R (t,t') = iort<t' 



df 2 



dt 



V ( ^{t,t',k) = 8{t-t') ; D A (t,t')=0 for t > t' 
V { °\t,t',k) = 



(7.37) 
(7.38) 
(7.39) 



Since a ~ T/a and in the initial state T > T c ~ n/y/X then for time scales in the intermediate regime [i 1 < t < t n i 
it is clear that t ^> o~ x . 



For k <C a and i > 1/<t (when we can neglect the second order time derivatives in ( 7.37 - 7.39 )) we find 



""" ''/■>- L-ll—1 6{t-t') , V { °\t,t',k) 



a 



8(t' - 1) 



V < g\t,t',k)=i 



L (t+t') 



V^(t,t',k) = e 



= «-£(*-*': 



/a. 



(7.40) 
(7.41) 



The homogeneous solution F(t, t', k) obeys the linear integral equation 



F(t,t',k)^V ( ^(t~t',k) + \ I dtiD { a } (t-t u 1e) I dt 2 n l {t 2 ) d(h - t 2 ) + U R (t 1 ,t 2 ;k) F(t 2 ,t' ,k) + t *-> t' 



which can be solved in perturbation theory. Writing F(t, t' , k) = F^(t, t' , k) + a F^(t, t' , k) + 0(a 2 ) we find 



(7.42) 



F ( -°\t,t',k)=V { °'(t-t',k) 



(7.43) 



F^(t,t',k)= / dhVg'it-t^k) dt 2 U l (t 2 ) S(t 1 -t 2 )+n R (t 1 ,t 2 ;k) V ( °'{t 2 ~t',k) + {t^t r ) (7.44) 



Obviously F^(t,t',k) gives the vacuum contribution to the magnetic field spectrum and must be subtracted. As 
discussed in the previous section both the tadpole and the retarded (bubble) self-energy are of order e 2 /A near the 
non-linear time, whereas n> is of order e 2 /A 2 . 

Hence one can show iterating the evolution equation in e 2 that the homogeneous contribution S R (t, k) is subleading 
by one power of A with respect to S B (t, k). 

Furthermore, since we are interested in the soft momenta regime because only modes in the spinodal band q 2 < 



H 2 <C rrij- increase, we may approximate 



T T 
l + n q ~n q ~ — ~ 

W q WIT 



l \q+k\ 



1 + n\g. 



W+k\ 



(7.45) 



The long- wavelength mode functions in the spinodally unstable band are approximately given by (5.15) for < t < t n i, 
hence in this approximation of long-wavelength and high initial temperature the magnetic spectrum 



S B (t,k) = e 2 [ q 2 (l - cos 2 : 

Jq 

(i+Ov+^i / dukVcitM^f^r^iti) 



t 2 

dh kV c {tM,k)f q {ti)f w+% Ati) - 



In 



+ n q {l + n w+% \) 



it, 



dHkVc^t^r^l^At!) 



n l n \q+k\ 



dtxkVcitM^r^f^ih) 
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simplifies to 

S B (t,k) = 96 °*J 2 e-^* J q 4 dqd(co S 9)(l-cos 2 d) J\^f q {ti)f w+ - H {ti)dt 



2 



(7.46) 

The integrals over momenta and angles can be done straightforwardly in the limit k -C fj,. By using the high 



temperature expression for the nonlinear time (5.19) we find that the integral is dominated by the upper limit and 
we obtain 



S B (k,t~t«)~ ^ /2Na /^ e-^ium (7 . 47) 

with the correlation length or domain size given by £(t) = but now with the nonlinear time given by t n i — 

m(^)/2/i. 

An important aspect of this result is that the factor e~ 2k t l" is cancelled by a similar factor in the time integrand, this 
is a consequence of the fact that at long time t > the integral is dominated by the upper limit. This cancellation 
in turn implies that the magnetic field spectrum generated by the instabilities is insensitive to the diffusion length 
Zdiff{t) ~ y/t/a- 

Besides some numerical factors and the arguments in the logarithms, the most important aspect as compared to the 
case of a quench in vacuum from the previous section is the factor Writing T = T C (T/T C ) and with T c = /i^/24/A 
we find 

a AT a T c \n(l/Na) 

'a ~ V24 NT « 1 ^ 

Therefore, the presence of a high conductivity plasma severely hinders the generation of magnetic fields. However, a 
noteworthy aspect is that up to the nonlinear time the magnetic field is still correlated over the size of the scalar field 
domains rather than the diffusion length Cdiff ~ \ftJo. This is an important point, the free field power spectrum for 
long-wavelengths k <C cr in a medium with high conductivity is given by 

S { °\t,k) = -i k 2 D%\t,t,k) = — e~^ (7.49) 

a 

which clearly displays the diffusion length scale = 2y/t/a. The diffusion length typically determines the spatial 
size of the region in which magnetic fields are correlated in the absence of non-equilibrium generation. The ratio 
between the domain size £(t) and the diffusion length scale £,diff(t) is given by 

^ - ^»1 (7.50) 



€diff(t) V A 4 



Where we have used the estimate (7.48). Thus an important conclusion of this study is that the magnetic fields 
generated via spinodal decomposition are correlated over regions comparable to the size of scalar field domains which 
are much larger than the diffusion scale. While this study does not directly apply to the cosmological situation, it is 
definitely encouraging and will be studied in more detail elsewhere [j39|. 

The spectrum for the electric field can be obtained from that of the magnetic field by simply replacing k T>c — > 2?c- 
In the soft regime and for time scales — <C t <C p- we have T>c — —k 2 /a 2 whereas kV c ~ kja. Therefore the electric 
field spectrum is suppressed by a factor k 2 /a 2 with respect to the magnetic field, namely 

S%{t,k) = £ s S%(t,k). (7.51) 

Thus in a high temperature plasma with large conductivity the non-equilibrium processes favor the generation of 
magnetic photons instead of electric photons, and again equipartition is not fulfilled. 

The energy density on wavelenghts A > L can be computed in closed form in the limits L ^> £(t n i) or L <§; £(t n i). 
The limit L ^ £(t n i) is the relevant one, if one is interested in the generation of large scale magnetic fields; the 
limit L <C £,(t n i) is important in order to estimate the power on small scales, this is important in cosmology since 
anisotropics in the cosmic microwave background imply severe constraints on the strength of magnetic fields at short 
scales. 
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Since we have used the expression for the mode functions in the spinodal band, small scales means scales much 
smaller than £ but still larger than jjT 1 . 



The ratio of the magnetic energy density on scales larger than L at the nonlinear time and the magnetic energy 
density in the radiation background, given by the Stefan-Boltzman law p 1 — 7r 2 T 4 /15 is given by 



While the prefactor may depend on the details of the cosmological setting, the factors (LT)~ 4 , (£T) -4 are purely 
dimensional and are ultimately the determinining factors for the strength of the generated magnetic fields on a given 
scale. These factors are invariant under the cosmological expansion and are determined by the ratio of the scales 
of interest today (galactic) to the thermal wavelength (today) of the cosmic microwave background radiation at the 
Wien peak. 

What can we learn for cosmology?: 

Although in this article we have focused on quenched phase transitions in Minkowski space-time, we expect that 
several of the results that we have found are robust and will emerge in a full and detailed study in cosmological space 
times. 

• The enhancement ~ 1/A 2 is expected on physical grounds since the current-current correlation function that 
determines the photon polarization is constructed of mode functions whose amplitude become non-perturbatively 
large near the end of the transition, at the nonlinear time scale. This is an unavoidable consequence of the 
classicalization of long- wavelength fluctuations during the phase transitional p8[ . 

• A quantity of relevance in the astrophysics of magnetic fields is the correlation length of the magnetic fields. 
A result that transpires from our analysis of the vacuum case and the case of a phase transition in a high 
temperature, or radiation dominated plasma, is that during the early and intermediate stages of the transition, 
the correlation length of the magnetic field tracks the size of the correlated scalar field domains. In the case of a 
highly conducting plasma we have seen that the size of domains is much larger than the diffusion length, hence 
the result obtained here is encouraging for generating magnetic fields on scales far larger than the diffusion 
length. 

• The large conductivity of the medium at high temperature hinders magnetogenesis. This result is in qualitative 
agreement with those in JlO| and implies that any calculational framework to obtain the spectrum of magnetic 
fields generated through non-equilibrium processes must account for the conductivity. In our study the conduc- 
tivity enters in the form /i 2 /<7 2 with \i the symmetry breaking scale, however in cosmology we expect the Planck 
scale to enter as well. 

• The important ratio r(L) that determines the viability of a dynamo to amplify a seed magnetic field features 
the factor (LT)~ 4 where L is the scale of interest. This factor is expected on dimensional grounds if the power 
spect rum at long wavelengths is cx k 2 as would be naively expected in a medium of high conductivity (see 
eq. ( |7.49| )). This factor is invariant under cosmological expansion and today LT ~ L/X cm ,b with A cm & is the 
wavelength at the Wien peak of the blackbody spectrum at temperature T, today T — T cm b. For a galactic 
scale (~ 30Kpc) LT ~ 10 24 and for the scale of galaxy clusters (~ IMpc) LT ~ 10 26 . Thus the ratio r(L) is 
naturally extremely small. Thus in order to overcome this enormous factor, which is obviously related to the 
cosmological expansion since the time of the phase transition, either the phase transition must occur fairly late, 
or the non-equilibrium processes must last for a long time or the generation of magnetic fields must involve a 
dramatic amplification factor. Of particular importance would be a phase transition in a radiation dominated 
cosmology that leads to a scaling solution for the dynamical evolution of the charged scalar fields. In this case 
the correlation function of the charged fields evolve in time in a scaling manner which will lead to a continuous 
generation of magnetic fields. This possibility is tantalizing because the 0{N) linear and non-linear sigma 
models do lead to a scaling solution in the large N limit |^9|, |37], 

We will report on our study of these possibilities in detail elsewhere [p9|. 



We find, 




(7.52) 




(7.53) 
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VIII. CONCLUSIONS 



In this article we have provided a framework to study consistently the generation of primordial magnetic fields 
produced by non-equilibrium phase transitions. The main premise is that during the process of spinodal decomposition 
which is the hallmark of non-equilibrium phase transitions (without free energy barriers) the instabilities that drive 
the process of phase separation and domain formation would lead to strong charge and current fluctuations if the 
scalar fields carry (hyper) charge. One of the main results of this article is a calculational framework derived from 
the underlying non-equilibrium Schwinger-Dyson equations. The result is an exact expression for the spectrum of 
magnetic and electric fields which is valid in all generality for scalar and spinor charged fields and in an arbitrary 



cosmological background [see eq.(6.45)]. The main ingredient is the transverse photon polarization and we have 
confirmed that this expression has the correct equilibrium limit. 

Separating the contribution from short wavelengths that are in local thermodynamic equilibrium from the long- 
wavelength fluctuations which fall out of equilibrium during the phase transition allows us to include the dissipative 
effects associated with the conductivity in a high temperature plasma. 

We studied the generation of magnetic fields during quenched phase transitions in Minkowski space time for a theory 
of N charged scalar fields coupled to an abelian gauge field (hypercharge) and one neutral scalar field (the inflaton). 
This is a reliable prelude to cosmology since the dynamics in Friedmann-Robertson- Walker backgrounds is similar to 
that in Minkowski space-time with a time dependent mass. Symmetry breaking occurs along the neutral direction 
and gauge symmetry is not spontaneously broken. We explicitly obtained the spe ctrum of ele ctric and magnetic fields 



to leading order in a (hypercharge) and to leading order in the large N [see eqs. (|7.23| )- (|7.28|) 1 . 

Two cases of cosmological relevance were studied, a transition in vacuum that models an inflationary phase tran- 
sition, and one from a high temperature phase to a low temperature phase that models a transition in a radiation 
dominated era. In the second case the high temperature plasma has a large conductivity. Separating the contribution 
from short wavelength fluctuations to the photon polarization we include the effects of the conductivity, while long 
wavelength fluctuations of the charged fields fall out of equilibrium during the transition and lead to the generation 
of magnetic and electric fields. We have provided explicit expressions for the magnetic and electric power spectra for 
long wavelengths [see eq.( 7.5I| )-(7.53)]. 



These studies have revealed several robust features which are expected to survive in cosmological spacetimes: 



The magnetic fields are correlated over length scales comparable to the size of the scalar field domains produced 
during the spinodal stage and in a high temperature plasma, this scale is much larger than the magnetic diffusion 
length. 

Equipartition between electric and magnetic fields does not hold out of equilibrium. For rolling down from 

2 

an initial vacuum state the ratio of electric to magnetic power behaves as ^ log A for long wavelengths [see 



eq. ( 7.28| )] . Electric photons overwhelm magnetic ones in this case. In the case of a conducting plasma, th e rati o 
between electric and magnetic power behaves as k 2 /a 2 for small k with a the DC conductivity [see eq.(7__5I)] 
therefore magnetic photons dominate over electric ones in this case. 

• The conductivity severely hinders the generation of long- wavelength magnetic fields. 

• The ratio r(L) between the energy density of magnetic fields on scales larger than L and the energy density 
of background radiation features an ubiquitous factor (LT)~ 4 which leads to large suppression factor which 
must be overcome in order for primordial seed fields to be amplified by the galactic dynamo mechanism. This 
point was originally made in ref . Q and is one of the formidable roadblocks that must be overcome towards the 
explanation of galactic and extragalactic magnetic fields from primordial seeds. 

We will report on our study of magnetogenesis in inflationary and radiation dominated cosmologies in a forthcoming 
article SI. 
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